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O. 
^ ■ Abstract 

Computational Group Theory is applied to indexed objects (tensors, spinors, 
and so on) with dummy indices. There are two groups to consider: one de- 
^T) ' scribes the intrinsic symmetries of the object and the other describes the 

interchange of names of dummy indices. The problem of finding canonical 
►^ , forms for indexed objects with dummy indices reduces to finding double coset 

^S| ' canonical representatives. Well known computational group algorithms are 

C^ . applied to index manipulation, which allow to address the simplification of 

1^^ I expressions with hundreds of indices going further to what is needed in prac- 

f^ ' tical applications. 
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o 
'^ • 1 Introduction 



Ref. |]1[ describes how Computational Group Theory provides tools for manipulating 
tensors with free indices. The tensors obey what we call permutation symmetries, 



Vh ' which are a set of tensor equations of the form 

where a{ii ■ ■ ■ z„) is a permutation of ii ■ ■ ■ in and e^ is either 1 or —1. This kind of 
symmetry can be described by finite group theory and the index manipulation can 
be performed using the algorithms of Computational Group Theory [|, ||, |, ||, ||, . 
The detailed description of symmetry as a group is given in ref. [Q. 

In this work we address the problem of applying Computational Group Theory 
for manipulating dummy indices. It is more complex then the free index problem, 
since one has to deal with two groups: the group that describes the symmetries 
of the indexed object and the group that describes the symmetry of interchange of 
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durniny indices. These groups act on a standard index configuration, generating sets 
of equivalent configurations. These sets are double cosets, which have already been 
studied in Computational Group Theory |^, |], ||. The most important concept 
for simplifying tensor expressions is the determination of canonical forms, which 
correspond to canonical representatives of single cosets for free indices and double 
cosets for dummy indices. The algorithms of the present work and the algorithms 
of refs. 0, flOl, and [|ll|| allow the manipulation of expressions built out of in- 
dexed objects obeying permutation symmetries, such as tensors, spinors, objects 
with gauge indices, and so on, with commutative or anticommutative properties. 
On the other hand, these algorithms do not solve yet the problem when there are 
algebraic constraints, such as the cyclic symmetry of the Riemann tensor. 

Manipulation of dummy indices can also be found in general algebraic expressions 
with sums and multiple integrals. For instance, the calculation of Feynman diagrams 
in Quantum Field Theory generates a large number of multiple integrals of the 
propagator which can in principle be reduced using the algorithms of this work by 
canonicalizing the integration variables. 

The structure of this paper is as follows. In section 2 we describe the representa- 
tion theory for dummy indices. In section 3 we describe the algorithm to canonicalize 
indexed objects with dummy indices, and in section 4 we discuss the algorithm com- 
plexity. In section 5 we discuss the simplification of general expressions in order to 
have a bird's-eye view of the problem, and present an example of canonicalizing a 
Riemann monomial of degree 3. 

We assume that the reader is familiar with the concepts and notations described 
in ref. [111. 



2 Representation theory for dummy indices 

Suppose that T is a fully contracted rank-2n tensor with symmetry S. We define 
the standard configuration as 

rpdi d2 d„ fr\\ 

■J- di d2 ' ' ' dn- y^) 

This configuration is associated with permutation +id, which is the least element of 
the symmetric group S'2„. Our first task is to determine the configurations that are 
equivalent to (|^). We know that the dummy index names can be interchanged. For 
example, the configuration 

T^'d/'d, ■■■'"d„ (3) 

is equivalent to (P) and is obtained by the action of the element +(1, 3)(2, 4) on @. 
This element is not in S in general. If the metric is symmetric, the configuration 

rp di d2 dn (A \ 



is also equivalent to (Q) and is obtained by the action of the element +(1,2). What 
is the group that describes these kinds of symmetries? Let D be a subgroup of 
H ® S2n generated by 

Kd = {+(1,2), +(3,4), ■■■ , +(2n- l,2n), 

+ (1, 3)(2, 4), +(3, 5)(4, 6), ■ ■ ■ , +{2n -3,2n- l){2n - 2, 2n)}, (5) 

with the base b^ = [1, 3, ■ ■ ■ , 2n — 1]. K^ is a strong generating set with respect 
to &£). The action of D on configuration @ yields all configurations that can be 
obtained from (0) by interchanging dummy index names or by using the symmetry 
of the metric. 

Besides the action oi D, we consider the action of S. If we take configuration 
(^ as the starting point, similar as we have done for the free index case, we have 
to apply D first, followed by S in order to obtain all configurations equivalent to 
(^. This order is crucial. If one applied an element of S first, the positions of the 
dummy indices would change and the application of D on this new configuration 
would make no sense. It would not be an interchange of dummy index names nor 
an interchange of contravariant index to a covariant one inside a pair. Let us see an 
example. Suppose that S = {+id, +(2, 3)}, and let us apply +(2, 3) on configuration 
(D followed by +(3,5)(4,6) G D. We obtain 

T"^ "' "' d, d, d, ■ ■ ■ "- d^, (6) 

which is not equivalent to (0) at all. The reverse order is perfectly fine, let us apply 
+ (3, 5)(4, 6) first, followed by +(2, 3): 

rpdi d-i d2 d„ /y^ 

-^ di ds d2 d„- \' ) 

The configuration above is equivalent to (0). 

The set (C) of all configurations equivalent to (0) is given by the action oi S x D 
on (D, i.e. 

C = {isd)iT'^ d/'d,--- '" dj, seS,deD}. (8) 

The set S X D is the double coset of S and D m. H ® S'2„, that contains the identity 
+id. The cardinality of this set is |S'||D|/|5' fl -D|. 

Consider a fully contracted configuration (Ti) that is not equivalent to (0), one 
can take (H) as an example. Suppose that Ti is obtained by acting ^f on (^). Then 
g ^ S X D. The set of all configurations equivalent to Ti is given by the action of 
the double coset S x g x D on @. The cardinality of this set is |S'||Z)|/|S'^ fl D|, 
where S^ is the conjugate set g~^ x S x g M. 



3 Algorithm to canonicalize tensors with dummy 
indices 

Now we address the following problem. Suppose one gives a fully contracted rank-2ra 
tensor which has some symmetry described by a set of tensor equations of the form 
(|l]). Find the canonical index configuration using the tensor symmetries, renaming 
of dummy indices, and the metric symmetries. 

In representation theory, this problem can be solved if one knows the solution of 
the following equivalent problem. Given a generating set Kg for the group S, and 
an element g = (e^,7r) E H ^ S2n, find the canonical representative of the double 
coset S X g X D, where D is the group generated by (|^) with respect to the base b^. 

Butler describes an algorithm for determining the double coset canonical 
representative for permutations groups. The input of his algorithm is: 

(a) a permutation group G acting on a set P with a base b = [61, ■ ■ ■ ,bk]; 

(b) subgroups A and B oi G given by a base and strong generating set; and 

(c) an element g of G. 

The algorithm determines the image of the base b under the first element of the 
double coset A x g x B. 

We modify Butler's algorithm in order to work within the direct product H®S2n- 
Fortunately, in the tensor problem, the subgroup D is fixed and we know beforehand 
a base and a strong generating set for it. Butler's algorithm keeps changing the base 
for the subgroup B (see item (b) above) during the determination of the image of 
the canonical representative. This base change is very simple for group D. 

Suppose that the settings in terms of tensor notation have already been converted 
into group notation. So, the input of the algorithm is: 

(a) n; 

(b) a base bs = [fei, ■ ■ ■ , h^] and strong generating set Ks for S; and 

(c) an element g = (evr, vr) G -ff S2n- 

The output is either the canonical representative g = (e^f , vr) of the double coset 
S X g X D or 0. The output occurs if and only if both (e^f, vf ) and (— e^, vf) are in 
S X g X D. To have the solution in terms of tensor notation when the output is not 
0, one simply acts ^ on (H). 

The algorithm basically consists of 2n — 1 loops. We describe the first two 
loops, which are enough to understand the whole process. The algorithm is formally 
described ahead. The base bs must be extended in order to be a base for S2n- Let 
bs = [bi, ■ ■ ■ ,bk, fefc+i, ■ ■ ■ , &2n-i] be the extended base and [pi, ■ ■ ■ ,P2n-i] the image 
of the canonical representative. The goal is to find pi,- ■ ■ ,P2n-i each at a time. 



First loop determines pi. 

The orbit b^^^^ gives all possible values of the first point of the image of the 
elements of the double coset S x g x D. Call this set IMAGESi. pi is the least point 
of IMAGESi with respect to 65- So, if Ab, = bf then 

IMAGESi = U i^. (9) 

The least point is calculated with respect to base 65. The order of points is 61 < 
■ ■ ■ < &2n- Before finding p2, we have to determine the pairs of elements {si,di) of 
(5*, D) that satisfy 

bl'"''=Pu (10) 

since, from now on, pi must remain as the first point. Suppose that (si, di) is a pair 
that satisfies ([TOD , then 

Ol =Pl- (11) 

So we have to determine a small set of pairs {si,di) and amplify it by using the 
stabilizers Sh^ and Dp^ in order to obtain all pairs {si,di) that satisfy (|T^). Notice 
that to determine Dp-^, we have to perform a base change so that pi becomes the 
first point of D . This is easily performed. The pairs {si,di) are stored in table 
TAB defined in the following way: 

TAB([2]) = (si,rfi), \i] e ALPHAi, (12) 

where ALPHAi is defined by 

ALPHAi = {[i],iG6fandz^Gj9f}. (13) 

The size of list [i] is given by the index 1 of ALPHAi. One can find ALPHAi using 

ALPHAi = Afc, n (pf)5"\ (14) 

We omit the dependence of {si,di) on the entries of ALPHAi, since the explicit 
notation (si([i]), (ii([z])) is cumbersome. It is important to keep in mind that for 
each entry of ALPHAi there is a correspondent pair (si, di). Note that the variables 
with index 1 are calculated in the first loop of the algorithm. 
The pair {si,di) corresponding to [i] is given by 

si = trace{i, us), 
di = trace{i^, i^d)~^, (15) 



where z/5 and ud are the Schreier vectors relative to the orbits of S and D respec- 
tively. First loop finishes here. 

For each pair {si,di) we calculate 62''^ ^ ^ ''Mn order to obtain 

IMAGES2= U ((&^)'^''^^)^''. (16) 

[i]gALPHAi 

IMAGES2 yields all images of 62 in the double coset S xgx D obeying the constraint 
(0). Then p2 is the least point of IMAGES2. 

Now we show how to find ALPHA2 and the associated pairs (52,(^2)- At this 
point, a pair (s2, ^2) has the following property: 

Kb2r''' = [pi,P2]- (17) 

For each pair (s2, ^2), define 

NEXT2 = (62'^)'^ n (pf"i)(^'^2)"\ (18) 

This is the set of images of 62 in S that yields p2 after applying g and d2. This set 
gives the points that extend ALPHAi. So 

ALPHA2 = {[i,j], [i] G ALPHAi, J e NEXT2}. (19) 

For each [i,j] in ALPHA2 we have to determine a pair (s2,c?2) that satisfies 

[61, 62]^"! •''2 ''''^'^'''^''1'^^ = [pi,P2]. (20) 

Let 

— 1 (2\ 

S2 = trace{f^ -.^s ) ^ ■^i' 
d2 = dix trace{f'^\v^^^)-\ (21) 

where Vg and v)^ are the Schreier vectors relative to the orbits of the stabilizers 
S*^^-* and D^'^^ respectively. We define the new entries of TAB as 

TAB([^,j]) = (S2,4), [^,j] e ALPHA2, (22) 

and clear the old ones. Second loop finishes here. 
In the ith loop, IMAGESj is given by 



IMAGES, = U {ihf'y^a'^^y , (23) 



where Si and di are obtained from TAB(L). NEXTj is given by 

NEXT, = {bf'y^ n (pf* 'V^*^"\ (24) 

and the pairs (sj, di) obey 

[bu■■■M^"'' = [Pu■■■,P^]■ (25) 

Now we present algorithm Canonical for dummy indices and the sub-routines Fi 
and F2. We use a pseudo-language that can be converted into programs of some 
computer algebra system. 

Algorithm Canonical (dummy indices) 
procedure double_coset_can_rep (n, bs, Ks, g) 

input: n number of pairs of dummy indices; 
bs = [^1, ■ ■ ■ ; bk] base of group S; 

Ks strong generating set of S with respect to base 65; and 
an element (7 = (e^r, vr) G if (X> S2n- 

output: g = (e^-, vf) canonical representative of the double coset S x g x D or 
if (e^, 7f) and (— e^, vr) are in S x g x D. 

begin 

bs '■= [&I5 ■ ■ ■ 5 bk, bk+i, ■ ■ ■ , &2ra-i] is the extension of bs in order 

to be a base for S2n', 
if metric is symmetric then 

Kd := {+(1,2), +(3,4), ■■■,+(2n-l,2n), 

+ (1,3)(2,4), +(3,5)(4,6), ■■■,+(2n-3,2n-l)(2n-2,2n)}; 
elseJf metric is antisymmetric then 

Kd := {-(1,2), -(3,4), ■■■,-(2n-l,2n), 

+ (1,3)(2,4), +(3,5)(4,6), ■■■,+(2n-3,2n-l)(2n-2,2n)}; 
else 

Ko := {+(1, 3)(2, 4), +(3, 5)(4, 6), ■ ■ ■ , +(2n - 3, 2n - l)(2n - 2, 2n)}; 
end if; 

6d:=[1, 3, ■■■ ,2n-l]; 
(* Initialize table TAB and ALPHA *) 
TAB{[]) ■.= {+id,+idy, 
ALPHA := {[]}; 

for i from 1 to 2 n — 1 do 

As :=all orbits of S (and calculate us - Schreier vector with respect to bs); 



A. := br^>- 

Ao := all orbits of Kj:,; 

(* IMAGES is given by eq. (^ *) 

IMAGES := map Fi on each entry of ALPHA passing TAB, 

Ad, Ab and g as extra arguments; 
Pi := least point of IMAGES with respect to base bs] 
6/5 := remove pj_i and move Pi (or pj — 1 if pi is even) to the 1st 

position in bf); 
ud := Schreier vector of Ko with respect to bo] 

for each L in ALPHA do 

s := 1st element of TAB(L); 
d :=2nd element of TAB(L); 
NEXT:= {AbY n A^f'^'^ ] 

for each j in NEXT do 

Si := trace{f , z/5) x s ; 
di := d X trace{j3'^,vo)~^ ', 
Li := append j to L; 
TAB(Li):=(si,c?i); 
end for; 
clear(TAB(L)); 
end for; 

ALPHA := indices of TAB that were assigned; 

(* Verify if there are 2 equal permutations of opposite sign in S x g x D *] 
if either Kg or Kd has some permutation with —1 then 
(* Calculate sgd for all (s, d) in TAB *) 
set_sgd:=map F2 on each entry of ALPHA passing TAB 

and g as extra arguments; 
if set_sgd has two equal permutations with opposite sign then 
break the loop and 
return 0; 
end if; 
end if; 

(* Find the stabilizers S^'+^^ and 0^^'+^^ *) 

Ks := remove permutations that have point hi from Ks'-, 

Kd := remove permutations that have point Pi from Kd] 



end for; 

g := F2(any entry of ALPHA); 
return g; 
end 



Sub routine Fi 

procedure Fi (L, TAB, Ad, Af,, g) 

input: L some entry of ALPHA; 

TAB table of elements (s, d); 
Af) all orbits of Kd; 
Ab orbit of some bi] and 

g = (e^, n) e H ^ S2n- 

output: ((Aft)*^'^)^'^^-^, where s and ci are the permutations associated with L. 
TAB(L) yields s and d. 

begin 

sgd ■= F2(L,TAB, (?); 

result := select the points of all partitions of A^, that 

have at least one point in (A^)*^'^; 
return result, 
end 

Sub routine F2 

procedure F2 (L, TAB, g) 

input: L some entry of ALPHA; 

TAB table of elements (s, d); 

g = (e^, n) e H ® S2n- 

output: sgd^ where s and d are the permutations associated with L. 

begin 

s := 1st element of TAB(L); 
d := 2nd element of TAB(L); 



result := s x g x d; 
return result] 
end 



4 Complexity 

The complexity of the general algorithm to find double coset canonical representa- 
tive is known to be exponential in the worst case 0, ^. On the other hand, the 
symmetries of tensor expressions are special cases of subgroups oi H ® S'2„, and ac- 
tual verifications show that in practical applications the algorithm is efficient. The 
symmetries of the Riemann tensor are one of the most complex that occur in prac- 
tice. Therefore, monomials built out of Riemann tensors are examples of complex 
tensor expressions. We have implemented algorithm Canonical and the auxiliary 
routines in Maple system [|l^ and have developed a program that generates at ran- 
dom Riemann monomials of any degree (number of Riemann tensors) with all indices 
contracted (Riemann scalar invariants). For each Riemann monomial we calculate 
the timing to find the canonical representative. We use a PC with a processor of 
600MHz. The vertical axis of the plot of Fig. 1 is the mean of 50 timings for each 
monomial. The horizontal axis is the degree. We have eliminated all timings of 
vanishing results. 



time(min) 




r O 15 

Riemann tensors 



Figure 1: Timing to find the canonical form of a Riemann monomial 
versus the degree. The dashed line is a fitting curve of the form y = 
9.3 X 10~^a;^, where x,y are the horizontal and vertical axis respectively. 
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/^Froin Fig. 1 one cannot prove that the algorithm is polynomial. It only shows that 
the implementation in Maple can handle monomials with large number of indices. 
The storage space is very low in order to produce the data. If we try to fit the 
experimental curve by a polynomial of the form y = ax^, for N < 5 the dashed 
curve passes above the experimental curve, and for A^ > 5 the dashed curve passes 
below for most of the points. The best polynomial using the least square method 
is y = 9.3 X 10"^ x^. Notice that the deviation from the polynomial curve depends 
on the degree due to the fact that 50 timings give worse and worse statistics with 
increasing degree. 

5 Canonicalization of General Expressions 

Consider an algebraic expression with indexed objects of tensorial nature. The 
product of these objects can be commutative or anticommutative. If we expand the 
expression, it becomes a sum of monomials. Refs. |T^ and Jill describe a method 



for merging monomials into single indexed objects, which inherit the symmetries of 
the original objects. The commutative or anticommutative properties are converted 
into permutation symmetries of the merged object. At the end, the problem of 
manipulating an expression reduces to the problem of dealing with single indexed 
objects with free and dummy indices obeying permutation symmetries. 

Without loss of generality, suppose that the merged object is a tensor T with p 
free indices and q pairs of dummy indices. We define the standard configuration as 

T*^-"^'^d,---''d,. (26) 

We do not distinguish contravariant free indices from covariant ones. If the original 
configuration has covariant free indices, we pretend that they are contravariant and 
proceed until the end, when the character of the covariant indices is restored. This 
means that there is no preference of putting contravariant free indices in front of 
covariant ones or vice-versa. The minimal order is dictated by the base of S which 
we do not know a priori, since it is built out by the strong generating set algorithm. 
This choice follows the criteria of least computational effort. 

All configurations of (^) taking into account sign changes are given by the ap- 
plication of elements in H®Spj^2q on (|26|). Suppose one gives an index configuration. 
The algorithms to canonicalize free and dummy indices can be applied in sequence 
on this configuration. The first step is the application of the algorithm of ref. |Q in 
order to find the canonical ordering and positions of the free indices. The next step 
is the application of the algorithm of section 3, translating the points [1, ■ ■ ■ , 2g] 
to the current positions of dummy indices. If [/i, ■ ■ ■ , hql are the new positions in 
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increasing order, then group D is strongly generated by 
Kd = { (^1, ^2)5 ■ ■ ■ 5 yhq-i^ hq), 

{h,h){h,U), ■■■ , {hq-3j2q-l){hq-2,hq)} (27) 

with respect to the base [h, h,- ■ ■ , hq-i], if the metric is symmetric. 

For example, let R"-'^'^'^ be the Riemann tensor and we want to canonicalize ex- 
pression 

Rd,d/"''Rds''"''Rd/'d/'. (28) 



Ref. ||T^ describes how this expression merges into a single tensor, which is 

Tb a d2 di d4 ds ds /'on^ 

ds d2 ds d4 di l-^yj 

with the following permutations symmetries 

Ks = {-(1,2), -(3,4), -(5,6), -(7,8), -(9,10), -(11,12), 
(1,3)(2,4), (5,7)(6,8), (9,11)(10,12), 

(5,9)(6,10)(7,11)(8,12)}. (30) 

Ks is a strong generating set. The standard configuration is 

T''"''d, ■■■"'d,. (31) 

The element oi H ® Su, which acts on the standard configuration (Rl|) and yields 



(71 = (1,12, 11, 4, 5, 6, 8, 9, 10, 7, 3). (32) 

Now we call the algorithm Canonical for free indices (ref. IH) with the following 
input: gi, Ks, and 65 = [1, ■ ■ ■ , 11]. We are using the simplest base in order to help 
the visualization of the order of the indices, and we are aware that it has unnecessary 
points. The output of the algorithm is 

^2 = -(2,5,6,8,9,10,7,3)(4,12,ll), (33) 

which corresponds to 

rpa d2 b di di di, ds ('iA\ 

-L ds d2 da di di ■ K'i^) 

The free indices are in the canonical positions, which are given by 

[1,2]^^-' = [1,3], (35) 

12 



and the positions of dummy indices are 

[3, ■■■,12] 5^ "' = [7, 11, 2, 5, 10, 6, 8, 9, 12, 4]. (36) 



Sorting with respect to the basis and concatenating (^), (p6D; converting to disjoint 
cycle notation we obtain 

/i = -(2,3), (37) 

which is the group element that converts K£, given by (^) to K£, given by (^^ via 
conjugation, i.e. K£, = h~^ x K^ x h. The input of the algorithm Canonical for 
dummy indices (section 3) is 

Ks,,, = {-(5, 6), -(7, 8), -(9, 10), -(11, 12), 

(5, 7)(6, 8), (9, 11)(10, 12), (5, 9)(6, 10)(7, 11)(8, 12)}, (38) 

65= [2, 4, 5, ■■■,11], (39) 

and 

93 = 92h = -(2, 5, 6, 8, 9, 10, 7)(4, 12, 11). (40) 

The algorithm must be modified so that the generating set for group D must be 

7^^ = 1(2,4), (5,6), (7,8), (9,10), (11,12), 

(2, 5)(4, 6), (5, 7)(6, 8), (7, 9)(8, 10), (9, 11)(10, 12)}, (41) 

with base bf) = [2, 5, 7, 9, 11]. The output is 

^4 = -(4,5)(6,7,9)(8,ll). (42) 

The permutations g^ and (74 do not act on the standard configuration (|3T|). They 
act on 

-^ di d2 ■ ■ ■ ds- l4dj 

The final answer is 

95 = 94h-' = -(2,3)(4,5)(6,7,9)(8,11). (44) 

In terms of tensor notation, the canonical form is 

-R-<iibd2j^^^d,d,d,j^^^^^^^^^^ (45) 

which is obtained acting g^ on (^) and splitting back the merged tensor. 
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